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A recently developed Standard-Model Extension (SME) formalism for neutrino oscillations that 
includes Lorentz and CPT violation is used to analyze the sidereal time variation of the neutrino 
event excess measured by the Liquid Scintillator Neutrino Detector (LSND) experiment. The LSND 
experiment, performed at Los Alamos National Laboratory, observed an excess, consistent with 
neutrino oscillations, of v £ in a beam of P M . It is determined that the LSND oscillation signal is 
consistent with no sidereal variation. However, there are several combinations of SME coefficients 
that describe the LSND data; both with and without sidereal variations. The scale of Lorentz and 
CPT violation extracted from the LSND data is of order 10~ 19 GeV for the SME coefficients ah 
and E x cl- This solution for Lorentz and CPT violating neutrino oscillations may be tested by 
other short baseline neutrino oscillation experiments, such as the MiniBooNE experiment. 



PACS numbers: 11.30.Cp, 14.60.Pq, 14.60.St 
I. INTRODUCTION 

Lorentz symmetry is one of the most fundamental ideas 
of both relativistic local quantum field theory and general 
relativity. Early tests, such as the Michelson-Morley and 
Kennedy-Thorndike experiments have established that 
Lorentz symmetry is an exact symmetry of nature. So it 
is natural to assume that Lorentz symmetry is an exact 
symmetry in the standard model (SM) of particle physics. 
However, since the SM does not address gravity, a funda- 
mental theory of Planck-scale physics (Mp ~ 10 19 GeV), 
including string theory pj and quantum gravity , may 
violate Lorentz and CPT symmetry 0. 

If limited to conventional relativistic quantum mechan- 
ics, it is possible to establish a self-consistent low-energy 
effective theory with Lorentz and CPT violation; this 
is called the standard-model Extension (SME) 4]. The 
minimal-SME formalism has all the conventional prop- 
erties of the standard model including observer Lorentz 
covariance, power counting renormalizability, energy mo- 
mentum conservation, quantized field, micro causality, 
and spin-statistics with particle Lorentz and CPT viola- 
tion due to background Lorentz tensor fields of the uni- 



* Deceased 



verse. The minimal SME also has SU(3) C x SU{2) L x 
U(1)y gauge invariance. Since the background Lorentz 
tensor fields are fixed in spacetime, by definition, they 
do not transform under an active transformation law. 
That implies rotation and boost dependence of physics 
in a specific frame. This formalism focuses on the inverse 
Planck-scale effect which is believed to be suppressed by 
at least one order of the inverse Planck mass (~ -M-, 
where E is the energy scale of the system under consid- 
eration). Therefore, the physics quantities involved in 
the formalism are perturbative. 

Surprisingly, atomic physics has achieved this sensitiv- 
ity level, and extensive experimental studies have been 
done (see, for example, Ref. 0). A recent experiment [I| 
of this type reaches a sensitivity to a specific combination 
of SME coefficients to order ~ 1CP 32 GeV, well beyond 
a basic estimate of the scale of new physics. In addi- 
tion, spectral polarimetry of distant cosmological sources 
yields a similar sensitivity for another combination of 
SME coefficients @. However, many of the SME co- 
efficients still have no experimental bounds. 

Similarly, quantum interference experiments, such as 
meson oscillations, are also sensitive to the small effect 
of Lorentz and CPT violation [7j . Tests have been made 
using data from many experiments, including KTeV 
FOCUS % BaBar [jjj, BELLE and OPAL 0. 
Recently, the SME formalism for neutrino oscillations, 
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another type of quantum interference experiment, has 
become available 
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Here, vacuum expectation values that contain Cj, B , d/^ B , 



II. THE LSND EVIDENCE FOR NEUTRINO 
OSCILLATIONS 



The Liquid Scintillator Neutrino Detector (LSND) ex- 
periment ^3 , completed at the Los Alamos National 
Laboratory (LANL), observed an excess of v e in a beam 
of Vn created from /i + decay at rest. The data analy- 

e + n events with 
If interpreted as 
to v e oscillations, this D e excess implies a two-neutrino 



sis used the sample of detected v e p —. 
positron energy 20 < E e + < 60 MeV. 



oscillation probability of (0.264 ± 0.067 ± 0.045)%. Here 
the first error is statistical and the second error is system- 
atic (neutrino flux, particle detection efficiency, cross sec- 
tions, etc.). Despite the evidence for neutrino oscillations 
from solar neutrinos [3 13 13 13 13 13, atmospheric 
neutrinos ^jl , accelerator neutrinos |23| , and reactor 
neutrinos |24|. the oscillation signal observed at LSND 
remains a puzzle. Since the neutrino sector is thought 
as likely to reveal new physics, the LSND anomaly is of- 
ten explained with new ideas such as a mass-difference 
CPT- violating model (see, for example Ref . |3] ) or ster- 
ile neutrino models (see Ref. 13] for a recent example). 
The MiniBooNE experiment [27J at Fermilab is currently 
taking data to test the LSND signal. 



III. LORENTZ VIOLATING NEUTRINO 
OSCILLATIONS 

Perhaps LSND is seeing the first signal of Planck-scale 
physics [28j. To describe neutrino oscillations, includ- 
ing Lorentz and CPT violation, a recently developed for- 
malism for neutrino oscillations 0, |3| using the SME 
framework |4j is employed. This framework allows for a 
sidereal time variation of the neutrino oscillation proba- 
bility. 

Within the SME framework, the neutrino free field La- 
grangian becomes, 
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The first term of and the first and second terms 
of Mab are the only nonzero terms in the case of con- 
ventional neutrino oscillations. The remaining terms in 
this Lagrangian represent the physics of the background 
fields. In general, the background Lorentz tensor fields 
are an infinite series, but if the focus is on a low-energy ef- 
fective theory, these eight additional fields are complete. 
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H^ab are CPT-even terms while 
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and b^g are CPT-odd by definition of the back- 
ground fields. Notice that each background field has fla- 
vor indices (A and B) that, unlike other systems, bring 
additional complication for the neutrino sector. 

This Lagrangian leads to the modified Dirac equation, 
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After some manipulation, this yields the effective 
Hamiltonian for active neutrino oscillations fl3| . In par- 
ticular, the effective Hamiltonian for active antincutrino 
to antincutrino oscillations is, 



(hcs)ab = \p\5 ab + {m 2 \l\ b 
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Here, the effective Hamiltonian is a 3 x 3 flavor Ma- 
jorana basis matrix of three active, right-han ded, an- 
tineutrinos. The original effective Hamiltonian [13j can 
describe v — u, v — v, and v — T> oscillations, but, in 
this work, lepton-number violating v — v oscillations are 
not considered. Therefore, the neutrino and antincutrino 
sectors can be diagonalized separately. There is some 
coupling of SME coefficients, (ol)^ = ( a )ab + (^)ab an< ^ 
(c L )^ = [cf a l + {df a v b . Also, other types of SME coeffi- 
cients do not show up in this analysis. For the usual con- 
ventional neutrino oscillation case, the effective Hamil- 
tonian (Eq. |SJ contains only the first two terms. Then, 
the neutrino oscillation probability depends on Am 2 and 
the mixing matrix. But, in this general form, including 
possible Lorentz and CPT violation, the diagonalization 
of the effective Hamiltonian is more complicated and, in 
general, it can not be represented by Am 2 and the mixing 
matrix alone. 



IV. THE SHORT-BASELINE APPROXIMATION 

If the baseline of the neutrino beam is short compared 
with the neutrino oscillation length, L, the neutrino os- 
cillation probability can be expanded with an effective 
Hamiltonian. Expressed to leading order in h e g |29| . 
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Note that in this equation, unlike the equations above, 
h and c have been explicitly included. Since, in the ef- 
fective Hamiltonian, contains information about the 
propagation direction of the neutrino, this oscillation 
probability depends on the neutrino propagation direc- 
tion. In order to form a phenomenological expression for 
the neutrino oscillation probability, it is most convenient 
to use a coordinate system fixed to the experiment [6ll30|. 
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FIG. 1: Coordinate systems for the sidereal time variation 
analysis: a) the Sun-centered system, b) the Earth-centered 
system, c) the LANL local coordinate system, and d) the 
definition of Tea = 0. 



The standard choice is a Sun-centered system (FigEk) 
that is, to a good approximation, an inertial frame for 
the experiment. 

With this choice of coordinates, the neutrino oscilla- 
tion probability becomes, 
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+ {A c ) eft cosw e T e + (B s )efi sin2w e T<j 
+(B c ) m cos2w e T e | 2 , 
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where cu& is the sidereal frequency (=27r/23h 56min 4.1s) 
and T@ is the sidereal time as measured from a standard 
origin. Note that Pa may depend on the sidereal 
time. 

These parameters, (C) S p, (A s )ep, (A c ) m , (B s )sp,, and 
(Bc)sp,y depend on the SME coefficients [a L Y and (c L y v 
and the neutrino propagation direction unit vectors N x , 
N Y , and N z in the Sun-centered system. The direction 
unit vectors depend on the colatitude \ °f the experiment 
in the Earth-centered system (Fig^) and the zenith and 
azimuthal angles 9 and cj> of the 9^ beam in the experi- 
ment local coordinate system (FigQJ;). 



cos x sin o cos <p + sin \ cos " 

sin 9 sin 4> I (8) 

- sin x sin 9 cos <j> + cos x cos C 



For neutrinos from the Los Alamos Neutron Science 
Center (LANSCE) beam to the LSND detector, x = 
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54.1°, 9 = 99.0° and = 82.6° |3J. The sidereal time is 
defined to be zero (T e = 0) at LANL local midnight on 
the autumnal equinox (FigGJi). At that time, the y axis 
of the Earth-centered system coincides with the Y axis 
of the Sun-centered system. The estimated error using 
this definition is three minutes, which is small compared 
to the time scale sensitivity of this analysis. 

Combining these values for neutrino propagation unit 
vectors with the detailed expression of t he p arameters, 
(C) g p,, (A s )ep,, (A c )efi, (B s )gp„ and (B c ) sfl :29J, yields, for 
the particular case of the LSND experiment: 
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rhgfi, has been in- 



In Eq. El the mass-squared term, 
eluded. This allows for conventional massive-neutrino 
oscillations in addition to the Lorentz-violation oscilla- 
tions. It is assumed that the size of this term does not 
invalidate the short-baseline approximation |29j . 



V. SIDEREAL TIME DISTRIBUTION OF THE 
LSND DATA 

In conventional explanations of neutrino oscillations, 
the oscillation probability is independent of sidereal time 
and, therefore, the sidereal time distribution of oscilla- 
tion events is expected to be constant. In the Lorentz 
and CPT violating model of neutrino oscillations consid- 
ered here, nonzero values of the model parameters could 
exhibit themselves as modulations to the sidereal time 
distribution (as in Eq. Q) . The sidereal time dependence 
of candidate oscillation events from the LSND data sam- 
ple has been examined and subjected to statistical tests 
to quantify any evidence for a sidereal variation. 

In the analysis of the final LSND data set Q, 205 
neutrino oscillation candidate events were reported with 
positron energy in the range 20 < E e + < 60 MeV and 
with an identified neutron-capture photon. There are two 
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classes of background in the oscillation sample: beam- 
unrelated and beam-related (^-induced). The beam- 
unrelated backgrounds arise from cosmic ray processes. 
It is measured in beam-off data and then subtracted from 
the beam-on data. The beam-related backgrounds are 
calculated from known neutrino (nonoscillation) interac- 
tions. 

The neutrino beam used for the LSND experiment 
was produced using protons from the LANSCE accel- 
erator |3^|. The proton beam was delivered at approx- 
imately 100Hz in pulses of 600 ^s duration. The de- 
tector was triggered independently of the state of the 
beam and the beam status was recorded. In this man- 
ner, beam-off data was taken continuously in the time be- 
tween beam pulses. The resulting beam-off data set was 
approximately 16 times larger than the beam-on data set. 
This allowed for an accurate measurement of the beam- 
unrelated background by weighting the beam-off data by 
the beam duty- factor (calculated for each run). 

The estimated number of beam-unrelated and v- 
induced background events in this sample are 106.8 ±2.5 
and 39.2 ± 3.1, respectively. These events were col- 
lected during experimental running in 1993 through 1998. 
There were six sets of runs, one in each of these years. 
The GPS (Global Positioning System) time stamp, nec- 
essary for this analysis, was not included into the LSND 
data stream until midway through the 1994 run period. 
Because of this, only 186 of the 205 oscillation candi- 
date events could be used in this analysis. The expected 
numbers of beam-off and ^-induced backgrounds in this 
smaller sample are 94.0 ±2.3 and 35.6 ±2.8, respectively. 

Ideally, an experiment to search for sidereal variations 
in a signal would run continuously throughout the calen- 
dar year so that one particular sidereal time bin would 
be drawn from the entire range of local time. This was 
not the case with LSND, but runs did cover the space of 
local time vs sidereal time with reasonable completeness, 
as can be seen in Figure [3 The Los Alamos (clock) time 
can be determined from Greenwich Mean (GM) time by 
subtracting 6 (7) hours in the summer (winter). 

To quantify the statistical significance of any side- 
real time variation in the data, we employed two dif- 
ferent statistical tests: a Pearson' s- y 2 test [33l Q an d 
a Kolmogorov-Smirnov (KS) test [34| • In both of these 
tests, the data were compared to the (null) hypothesis 
that the event rate is constant in sidereal time. Note 
that this null hypothesis is not that no oscillation sig- 
nal exists, but only that the signal is constant in time. 
We also examined the GM time distributions and applied 
these tests with the null hypothesis of an underlying dis- 
tribution that is constant in GM time. 

The Pearson's-x 2 (P-X 2 ); as implemented in this anal- 
ysis, is 
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where n is the total number of events in the sample, N is 
the number of time bins, and n; is the measured number 
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FIG. 2: Distribution of beam-on neutrino candidate events 
in a) run number vs. solar day (day 1= January 1) and b) 
GM time vs sidereal time. The year of each set of runs is 
indicated in a). 



of events in time bin i. The predicted number of events 
in each time bin, z/,-, is constant for each time bin. Note 
that this quantity is constructed with the variance of the 
expected number of events in the denominator. 

The P-x 2 statistic will follow, in the absence of sidereal 
time variations and with sufficient events per time bin, a 
X 2 distribution with number of degrees of freedom equal 
to the number of time bins minus one 0, 03] • The stan- 
dard criterion for sufficient events is that Vi > 5 [3^, . 
The binning for the beam-on data has been chosen to 
satisfy this. The p- value, (P(x 2 ), one minus the \ 2 cu- 
mulative distribution) can be extracted and interpreted 
as a confidence level that the null hypothesis explains the 
data. 

The KS test has the advantage that it works with un- 
binned data, thus eliminating the need to choose a bin- 
ning. It involves a comparison between the data and 
the null hypothesis via cumulative distributions. Unlike 
the P-x 2 test, it is sensitive to "runs" in the data, thus 
making the P-x 2 and KS tests complementary. The KS 
statistics reported here are the maximum cumulative de- 
viation, D n , and the KS probability, P(KS). The quan- 
tity P(KS), which is obtained from the known distribu- 
tion of D n jH, can be interpreted as a confidence level 
that the data is explained by the null hypothesis. 

In the LSND data set considered for this analysis, there 
were 1656 beam-off events passing the neutrino oscilla- 
tion cuts. These events, after weighting for the beam-on 
duty-factor, determine the number and distributions of 
beam-unrelated background events in the beam-on sam- 



5 



60 



1 



a 40 ± 



20 



60 



40 - 



20 



Till 



a) 



20000 40000 



60000 80000 
sidereal time (sees) 



20000 40000 60000 80000 

GM time (sees) 



t 

3 

c 



t 

3 

c 



10 



6 



2 - 



a) 



20000 40000 60000 80000 
sidereal time (sees) 




20000 



40000 



60000 80000 
GM time (sees) 



FIG. 3: Sidereal a) and GM b) time distributions of beam-off 
data using 37 time bins. The solid lines indicate the expected 
underlying distribution with no sidereal time variation. The 
errors bars displayed in this plot (and subsequent) are the 
square root of the number in each bin. 



FIG. 4: Sidereal a) and GM b) time distributions of the 186 
beam-on oscillation candidate events in 37 bins. The solid 
lines indicate the expected underlying distribution with no 
sidereal time variation. 



pie. The distribution of sidereal and GM times in 37 
time bins for these beam-off events is shown in Fig. 
The number of time bins chosen for this distribution was 
obtained by applying the N > 5 criterion for the beam-on 
data. The errors shown in Fig. (and subsequent figures) 
are the square root of the number of counts in the bin. 
Note that these errors are not used in the calculation of 
the P-x 2 (Eq.UUl. The P- X 2 is 29.6 for 37 sidereal time 
bins corresponding to P(x 2 ) — 0.77. The KS test on this 
same data yields D n = 0.019 and P(KS) = 0.60. These 
results indicate that the beam-off data are in reasonable 
agreement with the null hypothesis (no sidereal time de- 
pendence). The GM time distribution yields a slightly 
low P(KS) = 0.01, however, for this same distribution 
P(x 2 ) = 0.29. In addition, any GM time variations are 
distributed throughout a range in sidereal time. For these 
reasons, we conclude that there is no evidence for sub- 
stantial environmental or "day-night" sidereal variations 
in the beam-unrelated backgrounds. 

The sidereal and GM time distributions of the 186 os- 
cillation candidate events are shown in Figure 0] The 
P~X 2 for the sidereal time distributions is 44.8 for 37 
time bins. The corresponding p-value for the sidereal 
time distribution is P(x 2 ) = 0.15. A KS test applied to 
these distributions yields P(KS) = 0.234. The sidereal 
time distribution is slightly less compatible with no time 
variation as is evident in both of these statistical tests. 
However, the variation is not statistically significant. A 
KS test between beam-on and beam-off data was also 



applied and shows compatibility between the two data 
sets. The complete results from the statistical tests on 
the sidereal and GM time distributions for beam-on and 
beam-off data are summarized in Tabled 

To check the underlying assumption that the beam was 
delivered with equal efficiency throughout the sidereal 
day, a sample of 12 C(v e , e~) 12 N g , s , events was obtained 
by applying cuts to select for subsequent /9-decays of 
12 Ag. s . (as described in Ref. 35]). This procedure yields 
722 beam-on events with a beam-unrelated background 
of 17.5 events. The sidereal time distribution of these 
beam-on events are shown in Fig. [3] The P-x 2 for this 
sidereal time distributions is 29.3 for 37 time bins which 
corresponds to P(x 2 ) = 0.78. The Kolmogorov-Smirnov 
test yields D n = 0.020 and P(KS) = 0.94. The values 
indicate that the assumption of constant beam delivery 
averaged over the sidereal day and over all LSND runs, 
is consistent with the data. 



VI. THE EXTRACTION OF THE SME 
PARAMETERS 

While the LSND oscillation data examined in the pre- 
vious section shows no statistically significant sidereal 
time variation, it is interesting to examine the data in 
context of the SME model explained in Section Hvl First, 
this model does not require a sidereal time variation and, 
second, the LSND data set does allow for some sidereal 
time variation. 
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null hypothesis tests 




beam-on 
sidereal GM 


beam-off 
sidereal GM 


# of events 


186 


1656 


Pearson's \ 2 '- 


N bins 

x 2 
P(x 2 ) 


37 37 
44.8 27.6 
0.15 0.84 


37 37 
29.6 40.3 
0.77 0.29 


Kolmogorov-Smirnov: 


Dn 

P(KS) 


0.076 0.066 
0.234 0.386 


0.019 0.040 
0.604 0.010 



beam-on/beam-off tests 




sidereal GM 


D n 
P(KS) 


0.067 0.046 
0.432 0.864 



TABLE I: A summary of results from statistical tests on the 
sidereal and GM time distributions of the 20 < E e + < 60 MeV 
neutrino oscillation data. The null hypothesis tests com- 
pare the data with a constant time distribution. The beam- 
on/beam-off tests compare the two sets. 
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FIG. 5: Sidereal time distribution of beam-on data with 
12 C(u e , e~) 12 N g . B . cuts. The line indicates the average value. 

A maximum-likelihood method, with Eq. as a de- 
scription of the oscillation signal, was performed to 
extract allowed values of the SME parameters, (C)gp, 
(A)ejS) {A c )efi, (B s )sfi, and (B c )efi- In general, these 
parameters are complex — the special case is considered 
here where these parameters are real. Also, the values 
extracted are an effective average over the energy range 
of the LSND data set, 20 < E e + < 60 MeV. 

The parameters were extracted using an unbinned like- 
lihood function, 

-u N 



where N is the total number of events in the sample, fi s 
is the total predicted oscillation signal events, fib is the 
estimated number of background events, and fi — fi s +fib- 
The shape of the data in sidereal time is described with 
the functions J- s and T\>. T s depends on the SME pa- 
rameters as in Eq. [7|and Tb is assumed to be constant in 
sidereal time. The latter half of the likelihood function 
describes systematic errors on the signal and background 
events. In implementation, the natural log of the likeli- 
hood function, I (= In A), was used. Note that this func- 
tion describes both the shape and the overall number of 
events. 

Three different parameter combinations were consid- 
ered. 

• 1-parameter: 

(C)ep, 7^ 0) {•A s )efli (-^c)e/i; {Bsjejii ( B n )pn, = 

The "rotationally invariant" case [M 1371 l3Sl l39| . 

• 3-parameter: 

(C) m (A s ) m [AJsfi + 0; {B s )e fi , (B c ) m = 
Includes all of the CPT-odd terms of the minimal- 
SME model. 

• 5-parameter: 

(C)ep,i (-^-s)e/Zj (-^c)e/i5 {Bs*)efi: {Bc)efi 7^ 

Full minimal-SME model including both CPT-odd 
and CPT-even terms. 

Using each of these three parameter sets, the log likeli- 
hood, £, was calculated for the 186 candidate oscillation 
events as each of the parameters in the set was varied in 
a range around zero. The sidereal time for the parameter 
values that maximized £ is plotted together with the data 
in Figure El Note that the data in Figgis grouped into 
24 time bins instead of 37 as was used in Fig. 0] This is 
to allow for the quality of the fit to be more easily seen. 

The maximum-^ solutions are summarized below. The 
likelihood contours for the 3-parameter combination are 
shown in Fig.0 The (lcr) errors were calculated by deter- 
mining the parameter ranges where I decreased by 0.5 (1- 
parameter), 1.77 (3-parameters) , or 3.0 (5-parameters) 
from the maximum value. 

• 1-parameter: 

(C) m = 3.3 ±0.4 ±0.2 (16) 

• 3-parameter: 

There are two solutions within the lcr likelihood 
region (see Fig.[7J|. 
Solution 1 (maximum-^): 

(C) m = -0.2 ±1.0 ±0.3, 
(AJsp = 4.0 ±1.3 ±0.4, 
(Ao)sp = 1.9 ±1.8 ±0.4. (17) 
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FIG. 6: Sidereal time distribution of the LSND oscillation 
data in 24 time bins together with the maximum-^ solutions 
for the 1-parameter (solid line), 3-parameter (dotted), and 
5-parameter (dot-dashed) combinations. The dashed line in- 
dicates the estimated background contribution. 



Solution 2: 



3.3 ±0.5 ±0.3, 
0.1 ±0.6 ±0.2, 
-0.5 ±0.6 ±0.2. 



(18) 



• 5-parameter: Multiple (connected) solutions ex- 
ist in the 5-parameter case making a numerical ex- 
traction of errors impossible. The maximum- so- 
lution is: 





= -0.7, 


(-^s )e/i 


= 3.7, 


(-^c)e/i 


= 2.3, 


(B s )efi 


= 0.9, 


(Bc)ep 


= -0.6. 



(19) 



of 10~ 19 
are 



m 



GeV and 
the form 



All parameters have units 
the errors quoted above 
± (statistical) ± (systematic) . 

In all of these results, duplicate solutions exist with 
opposite signs for all of the parameters. Note that in the 
3-parameter case, the two solutions correspond to a large 
value for (C)gp with a small value for (A s )ep,, (A c )gp, and 
vice versa. The small- (C)g^, large- (A s ) g^, (A c )efi solu- 
tion is only slightly favored over the large-(C) g p, small- 
{A s )efi, {A c )efi solution. This is because the sinusoidal 
terms in Eq. improve the description of the data in 
sidereal time, although, an oscillation probability that 
is constant in sidereal time is consistent with the data 
(as was reported in Sec.EJ. Note also that the solution 
where (C)gp is the only nonzero term can be identified 
with the conventional neutrino oscillation description via 
the first term of Eq.0 A solution with all parameters « 
is highly disfavored. This is equivalent to the statement 
that the LSND oscillation excess is statistically signifi- 
cant. 
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FIG. 7: Log likelihood values for the 3-parameter description 
of the LSND sidereal time distribution: a) (As)sfi vs (C)ep., 
b) (Ac) e.;i vs (C)ep,, and c) (Ac) sp. vs (A s )ep.- The contours 
in a)-c) indicate the 1-a (total error) allowed regions (£ > 
-foiax — 1.77) and the stars indicate the maximum-^ parameter 
values. 



Since the oscillation probability depends on the SME 
parameters squared, the results for the SME parameters 
obtained above are more easily compared to the mea- 
sured oscillation probability from LSND via combina- 
tions of the squares of the parameters. The value re- 
sulting from the 1-parameter solution is 



1(C), 



10.7 ±2.6 ±1.3 (10" iy GeV) 



(20) 



The values for the parameter square sum resulting from 
the multiparameter combinations are more highly con- 
strained than for individual parameters. The value ex- 
tracted from the 3-parameter solution is 



\{C)^^\\(A s )- eri \ 2 + \\(A c U L \ 2 

= 9.9 ± 2.3 ± 1.4 (10~ 19 GeV) 2 , 

and from the 5-parameter solution, 

I (CM 



(21) 



2 I (-^ls)e/j. | + 2 l("^c)e/-i| 



+ i|(S s )g,| 2 ±i|(S c )g p | 2 

10.5 ± 2.4 ± 1.4 (10~ 19 GeV) 2 . 



(22) 



These results for the combination of SME parameters are 
consistent with the previously reported oscillation prob- 
ability from LSND [lj| and with the estimate presented 
in Ref. El. 
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FIG. 8: Sidereal a) and GM b) time distributions of the high- 
energy, beam-off data in 14 bins. The solid lines indicate 
the expected underlying distribution with no sidereal time 
variation. 



FIG. 9: Sidereal a) and GM b) time distributions of the 
73 LSND high-energy, beam-on oscillation candidate events 
in 14 bins. The solid lines indicate the expected underlying 
distribution with no sidereal time variation. 



VII. A HIGH-ENERGY SUBSET OF THE LSND 
DATA 

A high-energy subset of the LSND data with a positron 
energy cut, 36 < E e + < 60 MeV, is interesting to ex- 
amine separately. The i/-induced background is reduced 
in this sample Furthermore, the l/E prefactor to 

the mass term in Eq. would suppress the conventional 
oscillation terms relative to any Lorentz-violation terms 
present. 

This reduced data set consists of 73 beam-on events 
with expected beam-unrelated and ^-induced back- 
ground events of 31.3 ± 0.8 and 10.0 ± 0.8, respectively. 
The sidereal and GM time distributions of the 571 beam- 
off events passing these high-energy cuts are shown in 
Figure|HJ The P-x 2 is 12.2 for 14 sidereal time bins, cor- 
responding to P(x 2 ) — 0.51. The resulting p-value from 
the KS test to this distribution is P(KS) = 0.080. Again, 
these values show no reason to reject the null hypothesis 
for beam-off data. 

The sidereal time distribution of the high-energy 
beam-on data is shown in Figure |5J The P-% 2 is 20.4 
for 14 sidereal time bins corresponding to P(x 2 ) = 0.09. 
The resulting p-value for the KS test is P(KS) = 0.178. 
Although these values indicate a slightly reduced agree- 
ment with the null hypothesis, they do not indicate a 
statistically significant sidereal variation. The complete 
results from the statistical tests on the sidereal and GM 
time distributions for the high-energy data are summa- 
rized in Table ITT1 
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FIG. 10: Sidereal time distribution of the high-energy LSND 
oscillation data in 9 time bins together with the maximum- 
t solutions for the 1-parameter (solid line) and 3-parameter 
(dotted) combinations. The dashed line indicates the esti- 
mated background contribution. 



The maximum-likelihood procedure was applied to this 
high-energy data set using the 1- and 3-parameter com- 
binations described in Sect ion IVll The limited data sam- 
ple did not allow for the 5-parameter combination. Fig- 
ure llOl fwith a reduced bin size) shows 1- and 3-parameter 
maximum-! solutions superimposed on the high-energy 
data. Both parameter combinations produce acceptable 
descriptions of the data. 

The values for the parameter square sums extracted 
with the maximum-likelihood method are summarized 
below. The likelihood contours for the 3-parameter com- 
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null hypothesis tests 




beam-on 
sidereal GM 


beam-off 
sidereal GM 


# of events 


73 


571 


Pearson's \ 2 '- 


N bins 

x 2 
P(x 2 ) 


14 14 
20.4 9.7 
0.09 0.72 


14 14 
12.2 4.4 
0.51 0.99 


Kolmogorov-Smirnov: 


Dn 

P(KS) 


0.129 0.123 
0.178 0.221 


0.053 0.026 
0.080 0.826 



beam-on/beam-off tests 




sidereal GM 


D n 
P(KS) 


0.094 0.107 
0.621 0.451 



TABLE II: A summary of results from statistical tests on the 
sidereal and GM time distributions of the LSND 36 < E e + < 
60 MeV neutrino oscillation data. These values result from 
the same procedure as used for Table 



bination are shown in Fig. 1111 
• 1-parameter: 

|(C) gA | 2 = 10.7 ± 2.9 ± 1.5 (l(T 19 GeV) 2 (23) 



• 3-parameter: 

|(CM 2 + i|(AM 2 



2 |("4c)e/i| 



10.2 ± 2.7 ± 1.3 (10" 19 GeV) 2 



(24) 



As can be seen by comparing the results from the high- 
energy subset with the entire data set, there are no signif- 
icant differences. The time distributions from the high- 
energy subset are consistent with no sidereal variation 
and the results from the SME-parameter extraction are 
consistent with those obtained from the entire data set. 



VIII. A GLOBAL SOLUTION OF NEUTRINO 
OSCILLATIONS? 

To determine the implications of the allowed SME- 
parameter values extracted from the LSND data, con- 
sider the situation where the only nonzero term is (A s )ep,- 
This term was the largest in the maximum-^ solutions 
that allowed for sidereal variation. In this case, one or 
more of the SME coefficients, (oi)~, (oi)^, (ci)f^, 

( C i)g/Ti ( c L)sff an d { c l)X^ would be nonzero (as can be 
seen from Eq. I1U|I . 

A simple interpretation is that one of the a^-type SME 
coefficients is of order 10~ 19 GeV or one of the c^-type 
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FIG. 11: Log likelihood values for the 3-parameter descrip- 
tion of the high-energy LSND sidereal time distribution: a) 

(As)ep, VS (C)ep, b) (A c )gp VS (C)gp, and C) (A c )sp, VS (As)efi- 

The contours in a)-c) indicate the 1-a (total error) allowed re- 
gions (£ > fmax — 1.77) and the stars indicate the maximum-^ 
parameter values. 



is of order 10~ 17 (or E x cl ~ 10~ 19 GeV, where E 
is the neutrino energy). These values would have sig- 
nificant implications in other neutrino oscillation experi- 
ments and produce effects that have not been observed. 
In the simplest class of models, the acceptable maxi- 
mum scale of Lorentz and CPT violation for reactor neu- 
trino oscillations is ~ 10 -21 GeV, cl ~ 10~ 22 , and, 
for long-baseline neutrino oscillations, ~ 10 -22 GeV, 
c L - 10" 19 [H HI |H |H E3- However, there is no 
theoretical motivation that nature has chosen a simple 
solution for neutrino oscillations [40|. A global solution 
of neutrino oscillations with Lorentz and CPT violation 
that accommodates all the data may yet be obtainable 
within this SME framework. 

For this reason, it is important to search for sidereal 
variations in other short-baseline neutrino oscillation ex- 
periments. The data can be analyzed with the same 
method as presented here. 

The currently running MiniBooNE experiment |27| . 
with a different beam energy and with a beam, will 
be able to test these LSND solutions for Lorentz and 
CPT violating neutrino oscillations with a high-statistics 
appearance measurement. In particular, a measurement 
from MiniBooNE would provide an additional five con- 
straints (Eqs. I9I13|I on the SME coefficients. The neu- 
trino propagation vectors are different for MiniBooNE as 
the neutrino beamline is oriented toward compass north 
(as opposed to east for LSND). Also, the SME coeffi- 
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cients would be transformed for the neutrino case [29j- 
If MiniBooNE collects a significant set of data with a 
beam, an additional set of constraints with antineutrino 
coefficients would also be obtained. 

Of course, results from other neutrino oscillation ex- 
periments would add further valuable information. This 
has been investigated for Super-Kamiokande j^l) and MI- 
NOS 0. 

IX. CONCLUSIONS 

The neutrino oscillation candidate events from the 
LSND experiment have been examined for evidence of 
sidereal time variation — a possible signal for Lorentz 
violation in the neutrino sector. The oscillation excess is 
consistent with no sidereal time variation. An examina- 
tion of a high-energy subset of the data yields the same 
conclusion. 

A "smoking-gun" for Lorentz violation has not been 
found in the LSND signal. However, the data are ad- 



equately described within the SME neutrino oscillation 
formalism that includes both Lorentz and CPT viola- 
tion [T^L Efll ] . A maximum-likelihood method was used 
to determine allowed parameter regions for SME param- 
eter combinations. They indicate values on the order 
of 10~ 19 GeV for ar, and E x cl- These values are in 
the range expected for Planck-scale effects in the neu- 
trino sector. Future results from high-statistics oscilla- 
tion experiments will allow more stringent tests of the 
SME framework. 
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